Competition of Spin-Fluctuations and Phonons in Superconductivity of ZrZn 2 
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It has been long suspected that spin fluctuations in the weak itinerant ferromagnet ZrZn2 may 
lead to a triplet superconductivity in this material. Here we point out another possibility, a spatially 
inhomogeneous singlet superconducting state (a Fulde-Ferrell-Larkin-Ovchinnikov state). We report 
detailed electronic structure calculations, as well as calculations of the zone center phonons and their 
coupling with electrons. We find that the exchange splitting is nonuniform and may allow for gap 
formation at some parts of the Fermi surface. We also find that there is substantial coupling of Zr 
rattling modes with electrons, which can, in principle, provide the necessary pairing in the s-channel. 
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The recent discovery of superconductivity in the fer- 
romagnetic phase of ZrZn2 has revived interest in 
this unusual compound. ZrZn2 has traditionally been 
considered a prototypical example of a weak itinerant 
(Stoner) ferromagnet. Very small magnetic moments 
(values from 0.12 to 0.23 jib) have been reported. These 
do not saturate even at magnetic fields as high as 35 T, 
indicating softness of the magnetic moment amplitude 
and suggesting existence of soft longitudinal spin fluctu- 
ations. Already in the first report of ferromagnet ism 0] 
it was pointed out that ZrZn2 may be a ferromagnetic su- 
perconductor, of the type discussed earlier by Ginzburg 
H . Later the idea of triplet spin- fluctuation induced su- 
perconductivity in ferromagnetic ZrZn2 was theoretically 
elaborated by Fay and Appel ||]. However, experimental 
searches for superconductivity in ZrZn2 || were unsuc- 
cessful, till now. It is therefore very tempting to identify 
ZrZn2 as a triplet superconductor. This is supported by 
the fact that superconductivity seems to disappear with 
pressure at about the same point where the magnetic 
Curie temperature goes to zero. On the other hand, the 
low superconducting fraction does not allow accurate de- 
termination of the critical pressure for the superconduct- 
ing transition to clearly show that it coincides with that 
for the magnetism. Thus one cannot fully exclude the 
possibility that the two types of orders are spatially sep- 
arated. Still, it is likely that they do coexist, and if so, 
a key question is whether their interaction is construc- 
tive or destructive. A definite answer will require more 
study. However, before making any theoretical specula- 
tions one needs to get a detailed understanding of the 
electronic structure of this compound, and its relation 
with superconductivity and magnetism. Here we present 
an accurate analysis of the para- and ferromagnetic elec- 
tronic structure and identify what seems to be the most 
interesting and relevant features of the band structure for 
superconductivity. 

We find that, despite the very small magnetization, 
the rigid band model is inappropriate for ZrZn2. The 



exchange splitting, A xc (k), differs from band to band 
by more than a factor of two. Even more importantly, 
some bands have extremely small (Afc ~ 0.12 A -1 ) ex- 
change splitting in /c-space, making superconductivity in 
these bands less sensitive to magnetic ordering. Fur- 
thermore, the calculated Fermi surface shows substantial 
nesting features, suggesting that antiferromagnetic fluc- 
tuations may play a role in superconductivity. Finally, 
we found that the rattling Zr modes are soft and cou- 
ple strongly with electrons. We emphasize that, while 
triplet superconductivity is an interesting possibility, the 
facts are nevertheless compatible with an inhomogeneous 
Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state [Q of the 
s-wave symmetry as well. 




FIG. 1. Coordination of a Zr atom in ZrZn2. 



ZrZn2 occurs in a C15 crystal structure, where the Zr 
forms a diamond- type lattice, and the Zn forms a net- 
work of corner-sharing tetrahedra (another way to visu- 
alize this structure is to imagine a spinel without anions). 
An interesting aspect is the "over coordination" of Zr: it 
has 16 nearest neighbors consisting of 12 Zn, forming a 
truncated tetrahedron (Fig.|l|), at a distance (at T=50K 
H) 5.745 a.u., and 4 Zr at 6.001 a.u. The latter approxi- 
mately corresponds to the bond length in Zr metal. Since 
the metallic radius of Zn is 16% smaller than that of Zr, 
Zr and Zn do not form strong bonds. On the other hand, 
Zr, unlike carbon, does not form highly directional bonds, 
so 4 Zr-Zr bonds in ZrZn2 do not provide strong bonding 



either. This makes Zr "rattling" rather soft. At the same 
time, Zn has 8 neighbors at a distance 4.9 a.u., noticeably 
less than in Zn metal (5.04 - 5.51 a.u.). Correspondingly, 
one expect that Zn bond-stretching vibrations should be 
relatively hard. 




FIG. 2. LSDA band structure of ZrZn 2 . E F is at 0. 
Solid (dashed) lines show majority (minority) bands cor- 
responding to the magnetization of 0.2 /is/f.u. 

We calculated the band structure of both para- and 
ferromagnetic ZrZn2, using the experimental lattice pa- 
rameter of 13.858 a.u. |§. The calculations were done 
in the local spin density approximation (LSDA) || with 
the general potential linearized augmented planewave 
(LAPW) method. [Uj]. Well converged zone samplings 
and basis sets, including local orbitals Q] to treat the 
semicore states and relax linearization errors, were used. 
The spin-polarized band structure is shown in Fig.g. The 
Fermi surface is shown in Fig.J3|. Our calculations agree 
well with those published before (Refs. |13,|lj and refer- 
ences therein). The basic electronic characteristics of the 
nonmagnetic electronic structure are: density of states 
at the Fermi level, N(Ep) = 2.43 states/eV spin f.u., 
Fermi velocity, vf = 2.5 x 10 7 cm/sec, plasma frequency, 
uo p i = 4.0 eV, Hall concentration 5 x 10 22 e/cm 3 [ fJ4| . 
The bands near the Fermi level can be described in the 
first approximation as originating from the bonding com- 
bination of the t2 g Zr orbitals (the electronic sphere and 
the electronic rounded dodecahedron, both around the 
T point), and the bonding e g orbitals (the "intercon- 
nected pancakes"). The tubular network Fermi-surface is 
a hybridized combination of both types of states. How- 
ever, such a view is oversimplified. Comparing the width 
of the d-band of ZrZn2, and of the pure Zr sublattice 
(with Zn removed), we observe that the band width 
of the latter is about 30% smaller. Using the well- 
known nearest-neighbor tight binding formula, we see 

that W Zr Zn 2 ~ y&zr-zr + 12t Zr-zm which gives a 
rough estimate of tzr-Zn ~ 0.5tz r -Zr, n °t surprising, 
given the smaller radius of Zn. 

Although the effect of Zn on the total band width is 
relatively small (cf. with the fee Zr with the same Zr-Zr 
distance, which has twice larger bandwidth, Wzr-fcc ~ 

wl2t| r _ Zr flql), it is not neglegible, and, interestingly, 
it is nonuniform over the Brillouin zone. Specifically, the 



electron surface around the T point has more than 50% 
Zn character, while the other bands have nearly every- 
where less than 20% of it (Fig.|) |||. This leads, in turn, 
to nonuniform exchange splitting A^ c (k)(Fig.g|). For the 
Fermi surface near the T point, the splitting of the Fermi 
surface is extremely small, SkF ~ 0.017 A -1 . The FFLO 
period, 27r/5kF ~ 360 A, is long enough to allow for 
superconductivity in this band with a coherence length 
f ~ 290 A (according to Ref. 0). This band has an 
exceptionally small 5kp, but the other bands have small 
8k p as well. In fact, for the most of the Fermi surface 
5kp ^ 0.07 A -1 , corresponding to 2ir/5kF ~ 90 A. 

The spin-polarized bands, used to produce the lower 
panel of Fig.[|, are shown in FigJ|. These were calcu- 
lated by fixing the total magnetization to be 0.2 /ig/f.u. 
Fully relaxed calculations converge to the value of 0.7 
\±b /f.u., nearly 4 times larger than the accepted experi- 
mental value of 0.17-0.20 /i#/f.u. This is somewhat un- 
expected, for usually correlation effects, not completely 
accounted for in LSDA, tend to increase the tendency to 
magnetism. It may be that the actual samples are spa- 
tially inhomogeneous, with magnetic and nonmagnetic 
(and possibly superconducting) regions. Alternatively, 
the LSDA may overestimate tendency to magnetism in 
this material for some unknown reason. 

In vicinity of an itinerant ferromagnetic critical point, 
theory predicts a triplet superconductivity, with T c first 
growing as one moves in either direction from the critical 
point, and then decaying as the spin fluctuations become 
weaker in either nonmagnetic, or ferromagnetic phases. 
Such a behavior has been observed on the ferromagnetic 
side, for instance, in UGe2, and it was suggested that 
the difference between the longitudinal spin fluctuation 
spectra in the para- and ferromagnetic phases makes su- 
perconductivity much weaker and thus unobservable on 
the paramagnetic side |17]]. So, a triplet state remains a 
plausible explanation for the superconductivity in ZrZn2. 
We want to call attention, however, to another possibil- 
ity, namely an FFLO inhomogeneous s-wave supercon- 
ducting state. Both models can successfully explain the 
major experimental observations [[jj: high sensitivity to 
sample purity (impurity scattering would make the su- 
perconducting electrons from a band with a small ex- 
change splitting feel a large splitting of the other bands), 
absence of a specific heat jump (the Fermi surface sec- 
tions with a large splitting may have vanishing gaps), 
disappearence of superconductivity near the FM critical 
point (because of the pair-breaking effect of spin fluctua- 
tions). An apparent difference between the two models is 
that in the latter superconductivity not only is restored 
at pressures higher than critical (this is true for most 
models of the triplet superconductivity as well), but also, 
generically, T c grows much faster on the paramagnetic 
side than on the ferromagnetic one. The high pressure 
data reported in Ref. M do not extend far enough be- 
yond the critical point and do not give a definite answer 



to whether superconductivity reappears at higher pres- 
sures or not. Unless the strongly interacting Zr rattling 
modes would harden substantially, and/or lose their cou- 
pling with electrons, under pressure, the FFLO scenario 
suggests superconductivity re-entrance at a pressure high 
enough to suppress the spin fluctuations. We did check 
the pressure dependence of the T2 g mode, and found that 
under a 3% volume compression it stiffens by 11%, while 
it softens by 6% under a 3% expansion. 

A major question that arises in connection with this 
alternative is, what would be the pairing interaction be- 
hind the assumed s-wave state? It needs to be sufficiently 
strong to make an FFLO state energetically competitive. 
To gain more insight in that, we performed lattice dy- 
namics calculations for the zone center optical modes 
[pq| . There are 15 such modes at 6 distinctive frequen- 
cies: 4 triple-degenerate modes of T\ u (two), T2 n , and 
T^ g modes, one double-degenerate E u mode, and one 
non-degenerate A\ u mode (Table 1). The two softest 
modes are particularly interesting. Physically, they cor- 
respond to "rattling" motion of Zr inside the Zn cage, 
hence their softness and the small difference in their fre- 
quencies. Since these two modes join at the zone bound- 
ary, one can see that they have little dispersion, as is 
common for rattling modes. The even (T2g) mode is the 
only one that can couple with phonons at the Y point; 
but as these two modes are essentially local vibrations, 
they presumably couple with about the same strength 
for general wave vectors, so the coupling constant com- 
puted for a T2 9 mode should be a reasonable estimate 
for the coupling constants of all six rattling T2 9 and T\ u 
phonons. 
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FIG. 3. Fermi surface of ZrZn 2 . The three upper figures 
show the nongnatic Fermi surfaces, colored by the relative 
percentage of Zn character. The color bar goes from 5% 
to 60%. The bands shown are, from up down: bands 1 
and 2 (r-points at the corners), bands 1 and 3 (L-points 
at the corners), bands 1 and 4 (r-points at the corners). 
The bottom panel shows the Fermi surface 100 (k x — 0) 
cross-section. The widths of the lines correspond to actual 
exchange splitting for magnetization M = 0.2 \±b. 

The coupling constant, defined as X ep = 2 Y2k a ^( £ k« — 
E F )(dck a /du) 2 /2Mu 2 , where ^ ka S(e\^ a — E F ) is the 
density of states at the Fermi level per spin, M is the 
ionic mass, uo is the phonon frequency, and (de\^ a / du) 2 , 



is X, 



0.115. Note that this mode is triple degen- 



erate at T, so the total contribution to X ep from this 
mode is ^0.35. However, since we expect that all rat- 
tling modes of Zr would couple at a general point in 
the Brillouin zone with a similar strength, we estimate 
the total contribution from the six rattling modes to be 
~ 0.7. Finally, the remaining 12 modes (especially the 
two rather soft E u modes, which we expect to couple less 
with electons, should bring a smaller, but finite contribu- 
tion to the the total A ep , so we take as a rough estimate 
X ep ~ 1. With a Coulomb pseudopotential /i* = 0.1, 
the McMillan formula yields T c « 10 K. This is a rea- 
sonable number since spin fluctuations should reduce T c 
(and eventually eliminate superconductivity) in this sce- 
nario. This is also consistent with the transport coupling 
constant extracted from the linear part of the resistivity 
vs. temperature dependence [[L9[ of ZrZn2 single crystals 
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1.4. Assuming that A^, 



Ae» + A es , where 



A es characterizes coupling with spin fluctuations, we get 
A es ~ 0.4, and the modified McMillan formula then gives 
T c = (u; ph /1.2)exp[-1.02(l+X ep +X es )/(X ep -X es -jl*) ~ 



0.3 K. Here/i* = /x*(l + 0.62A ep + 0.62A es ). On the other 
hand, the specific heat coefficient of 47 mJ mol _1 K -2 , 
reported in Ref. Jy , corresponds to mass renormalization 
4.1, considerably larger than (1 + X ep + A es ) = 2.4. |pl| 

Another interesting observation is that the topology of 
the Fermi surface changes so sharply with the chemical 
potential, that the calculated coupling constant X ep is 
very sensitive to the exact position of the Fermi level in 
the band structure. In Fig.|| we show its dependence on 
the Fermi level position. We observe that if the Fermi 
level were 20 meV, X ep would be 50% larger. This unsual 
sensitivity is only partially accounted for by the structure 
of the density of states (Figfl). Of course, the tendency 
to magnetism also strongly depends on the position of 
the Fermi level. 




FIG. 4. Electron-phonon coupling constant, A ep , for 



one T2 9 mode as a function of the Fermi level position. 
Solid line is A ep , dashed line is density of states. 

In conclusion, we report accurate calculations of the 
electronic structure and the zone-center phonons in 
ZrZn2 both in para- and ferromagnetic state. The 
recently observed superconductivy is compatible with 
any of the three scenarios: (1) sample inhomogene- 
ity leading to separation of superconductivity and mag- 
netism in real space, (2) coexistence of ferromagnetism 
and inhomogeneous s-wave superconducting state (Fulde- 
Ferrell-Larkin-Ovchinnikov state), and (3) triplet spin- 
fluctuation induced superconductivity. Further exper- 
iments should be able distinguish between these three 
scenarios, since they predict rather different thermody- 
namic behavior: an exponential BCS-like temperature 
dependencies vs. a finite residual density of states vs. a 
power law due to the gap nodes. 
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Zone center phonon modes in ZrZn2 
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Zr 

120 
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133 142 182 250 


Zn 
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a Out-of-phase Zr rattling; b In-phase Zr rattling; c Zn 
breathing (the tetrahedra breathe out of phase with each 



other); d Mixed mode; e Zn breathing (all three tetrahedra breathe in phase). 



